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Abstract 

Owing to the demand for sustainable sex-control protocols in aquaculture, research In tilapia sex determination is gaining 
momentum. The mutual influence of environmental and genetic factors hampers disentangling the complex sex 
determination mechanism in Nile tilapia {Oreochromis niloticus). Previous linkage analyses have demonstrated quantitative 
trait loci for the phenotypic sex on linkage groups 1, 3, and 23. Quantitative trait loci for temperature-dependent sex 
reversal similarly reside on linkage group 23. The anti-Mullerian hormone gene {amh), located in this genomic region, is 
important for sexual fate in higher vertebrates, and shows sexually dimorphic expression in Nile tilapia. Therefore this study 
aimed at detecting allelic variants and marker-sex associations in the amh gene. Sequencing identified six allelic variants. A 
significant effect on the phenotypic sex for SNP ss831884014 (p<0.0017) was found by stepwise logistic regression. The 
remaining variants were not significantly associated. Functional annotation of SNP ss831884014 revealed a non-synonymous 
amino acid substitution in the amh protein. Consequently, a fluorescence resonance energy transfer (FRET) based 
genotyping assay was developed and validated with a representative sample of fish. A logistic linear model confirmed a 
highly significant effect of the treatment and genotype on the phenotypic sex, but not for the interaction term (treatment: 
p<0.0001; genotype: p<0.0025). An additive genetic model proved a linear allele substitution effect of 12% in individuals 
from controls and groups treated at high temperature, respectively. Moreover, the effect of the genotype on the male 
proportion was significantly higher in groups treated at high temperature, giving 31% more males on average of the three 
genotypes. In addition, the groups treated at high temperature showed a positive dominance deviation (-1-1 1.4% males). In 
summary, marker-assisted selection for amh variant ss831884014 seems to be highly beneficial to increase the male 
proportion in Nile tilapia, especially when applying temperature-induced sex reversal. 
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Introduction 

Tilapias, famed for their rapid and efTicient growth, their low 
position on the food chain and deemed the "aquatic chicken", 
provide a possibility to nourish the poor and to conquer export 
markets. Even though intensive aquaculture farming resulted in 
rapid breeding and stunted populations, hybridisation or applica- 
tion of synthetic hormones (17a-methyl testosterone) allowed the 
production of mono-sex male broods [1]. The latter technique can 
effectively prevent fry production before harvest. Showcasing 
tilapia as the first aquaculture fish species certified by the 
Aquaculture Stewardship Council [2], this emphasizes the need 
for more sustainable sex control protocols. As such, a temperature 
treatment of the fry might be a future alternative [3], [4]. The 
complex sex determination (SD) in Nile tilapia is comprised of an 
interaction between genetic (major and minor factors) and 
environment-Ztemperature-dependent factors [5]. Due to its 



complex SD, Nile tilapia is a well-suited model species for the 
evolution of genetic or environmental sex-determining mecha- 
nisms. Although an increasing number of genomic tools has been 
developed for Nile tilapia, such as a genetic map with formerly 23 
linkage groups (LG) [6], which has recently been resolved in 22 
LGs in a high-density radiation hybrid map with 1358 genetic 
markers [7] , BAG end sequences [8] , and an assembled draft of its 
genome (http://www.broadinstitute.org/ scientific-community/ 
science/ projects/ mammals-models/ vertebrates-invertebrates/ tilapia- 
/tUapia-genom), a comprehensive understanding of the SD 
mechanism is stiU outstanding. On the one hand, linkage analysis 
demonstrated some success, with QTL found for phenotypic sex on 
LG 1, 3, and 23 in intraspecific and interspecific crosses of tilapia 
species [9-1 7] . On the other hand, neither assigning vertebrate SD 
candidate genes to tilapia linkage groups [18-20] nor the analysis of 
expressed sequence tags for ovary- and testis-specific libraries 
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resulted in appreciable success so far [21]. Beside major QTL on 
LGl and LG3 [22], Eshel et al. 2012 recently fine-mapped a SD 
region on LG23 [15]. The QTL mapped to 13^0 cM and more 
precisely peaked at 22 cM (F=78.7; P<7.6xlO"'*). Linkage 
mapping located the SD QTL between markers GM597 and 
AR0124. The QTL region, located in scaffold 101 of the Nile 
tilapia genome assembly, harbours the anti-MiiUerian hormone 
gene (amh), among others. A further study proposed that 
temperature-dependent sex reversal similarly seems to be influenced 
in at least some families by an even broader QTL region (between 
GM283 and UNH898) on LG23 [23]. 

The amh gene is involved in the development of the urogenital 
system during embryogenesis by suppressing the Miillerian ducts 
of mammals, birds and reptiles [24—26]. Both amh and its target in 
the TGF'P pathway, the type II amh receptor (amhrll), are 
expressed in Sertoli cells of teleost [27]. Moreover, amh has been 
cloned in a number of fish species despite the fact that they lack 
Mullcrian ducts [28-31]. The expression of the amh gene can be 
detected starting at 3 days post fertilisation (dpi) in XX and XY 
gonads [15]. The authors also demonstrated a significantly higher 
expression of amh in all-male vs. all-female batches at 3 dpf 
Differences became more pronounced until 7 dpf [15]. Earlier 
investigations tended to show dimorphic expression in the gonads 
at later stages, i.e. after dmrtl at 19 dpf [32]. Moreover, sustained 
up regulation of dmrtl and amh during temperature sex-reversal 
from 13-15 dpf onwards and down regulation offoxl2 and/or 
cypl9ala at 17-19 dpf were observed during sexual development 
of Nile tilapia [33]. Even in the brain, amh expression appeared to 
be clearly sexually dimorphic during the period from 10 to 15 dpf 
[34]. 

This study reports the detection of allelic variants in the coding 
and 5' upstream region of amh. The development of a high- 
throughput genot^ping assay for a single SNP and its association 
with the phenotypic sex let amh appear as a candidate for 
autosomal and temperature-induced sexual development in Nile 
tilapia. 

Materials and Methods 

Stocks and Cross Design 

The aim of the present study was to identify alleles in the amh 
gene with an influence on autosomal and temperature-dependent 
sex reversal. For this purpose genetically all-female (XX) 
populations were used, which had been tested for their chromo- 
somal sex (XX vs. XY) via progeny testing (Figure 1). The 
genetically all-female population was derived from crosses between 
females from a selected line for low responsiveness to temperature 
(<60% males after 36°C treatment from 10-20 dpf) and 
temperature-sex reversed males from a line for high tempera- 
ture-responsiveness (>90°'!) males after 36°C treatment from 10- 
20 dpf) [3,4]. The genetically female line was developed through 
mating of 9 scx-rcvcrsed (AXX) sires derived from 7 high-line 
families, which were treated at high temperature, to 7 females 
from 6 low-hue families in order to produce a total of twelve 
families. One line cross family was derived from a within line 
mating (sire 105, dam 32, Table SI in File SI). Briefly, progenies 
were obtained through artificial reproduction (Figure 1). Fertilised 
eggs were incubated for 10 days at 28°C. After yolk sac absorption, 
larvae from each family were randomly distributed into two 
groups (n~110 larvae each). Temperature in the control groups 
was 28°C throughout the experiment, whereas treatment groups 
were kept at 36±0.5°C from 10-20 dpf according to [3]. From 20 
dpf onwards, control and groups treated at high temperature were 
raised in separate tanks at 28°C for at least 2 months. From each 



family treated at high temperature 15 adult males and 15 adult 
females were phenotyped alive, assessing the urogenital papilla 
and the types of gametes each fish produced. AH individuals that 
could not be confidently phenotyped, even after repeated 
inspection of the genital papilla, were excluded from further 
analysis. The remaining fish from the groups treated at high 
temperature and control groups were subjected to a lethal dose of 
phenoxyethanol (ethylene glycol monophenyl ether at 500 nL/L) 
and were immediately exanguated and dissected. The latter fish 
were phenotyped for their sex based on microscopical inspection 
of squashed gonads according to [35], classifying them into either 
testes or ovaries. 

Tissue samples from all individuals were collected from caudal 
fins and stored at — 20°C until DNA extraction. The DNA was 
isolated from a sample of finclips by phenol-chloroform extraction 
[36]. The sample was comprised of 305 individuals from 7 out 12 
line cross families including all 9 AXX founder sires of the 1 2 line 
cross families, which were all treated at high temperature. 
Furthermore DNA was extracted from 125 control individuals 
from 10 out of 12 line cross families including the 7 founder dams 
of the 12 line cross families and two additional females with no 
offspring. The numbers of fish genot^ped from each family and 
corresponding control or treatment group are gi\'en in Table SI 
and S5 in File SI. Water parameters measured during the 
experimental period were within the following range: oxygen > 
5.5 mg/L; pH 6.5-7.5; NH*+<0.5 mg/L; NO2"<0.25 mg/L. 
For first feeding, the fish were provided a diet rich in protein three 
times a day ad libitum (Tetra Werke, Germany; crude protein 
48.5%). Generally, from day 20 until day 60, fish were fed ad 
Ubitum with a trout feed (Skr(;tting FO Aqua Brut, Norway; crude 
protein 48.5%) while from day 60 the fish were provided a carp 
feed (Skretting C2 Pro Aqua K18, Norway; crude protein 36%). 
All procedures were in strict accordance with the recommenda- 
tions in the Guide for the Care and Use of Laboratory Animals of 
the German Animal Welfare Act [37]. This study was approved by 
the Institutional Animal Care and Use Committee of Goettingen 
University. 

Sequencing of the amh gene 

The genomic sequence of the amh gene was derived from 
Scaffold GL83 1234.1 of the Nile tilapia genome sequence 
deposited in the Ensembl database (http://www.ensembl.org; 
Orenill.O GCA_000 18823 5.1; location of amh: Scaffold 
GL831234.1, 1.688.687-1.691.779). Gene-specific primers were 
designed to cover the coding sequence (cds), including intron and 
5'UTR regions of the amh gene using the Primer3 software (see 
Table S2 in File SI). Each forward or reverse primer for 
sequencing of the amh gene was tailed at the 5'end with the 
Ml 3 universal forward or reverse primer to enable direct 
bidirectional sequencing on an 3130xL Genetic Analyzer (Applied 
Biosystems, Germany) using the BigDye Terminator v3.1 Cycle 
Sequencing Kit (Applied Biosystems). DNA for sequencing of the 
amh gene was derived from 93 individuals comprised three 
families including the corresponding three dams and sires (Table 
S6 in File SI). PGR was carried out using 20 ng of genomic DNA, 
Ix PGR buffer containing MgCl2, Ix Q;Solution, 10 pmol of 
each primer, 10 mM dNTPs and 2 U FastStart Taq DNA 
polymerase, in a final volume of 25 |J,1. All PCR components 
except the primers (MWG, Germany) and the 1 x Qjsolution 
(Qiagen, Germany) were purchased from Roche Diagnostics 
(Mannheim, Germany). PCR was performed using a Biometra T- 
3000 Thermocycler (Goettingen, Germany) with an initial 
denaturation at 95°C for 10 min, followed by 35 cycles of 92°C 
for 30 s; 60°C for 30 s and 72°C for 1 min with a final extension 
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Figure 1. Experimental design to obtain a genetically female (XX) population in order to study temperature effects on the male 
proportion in a selected line of Nile tilapia. 

doi:10.1371/journal.pone.0104795.g001 



at 72°C for 5 min. The fragment identity was controlled via gel- 
electrophoresis on 1.5-2% agarose gels. PGR products were then 
purified with Exo-SAP-IT (USB, Germany). The obtained 
sequences were trimmed, contigs were built, and SNPs were 
manually identified using the program software suite DNASTAR 
LasergeneG (DNASTAR, Inc., Germany). 

Genotyping of SNP ss831884014 within amh 

Gomparative sequencing revealed six SNPs within the amh gene 
(Figure 2). To genotype SNP ss831884014, exhibiting the largest 
effect on the phenotypic sex, a fluorescence resonance energy 
transfer (FRET) was developed using a LightGycler 480 instru- 
ment (Roche Diagnostics, Mannheim, Germany). Primers Fret-for 
and Fret-rev were designed to flank SNP ss831884014 yielding a 
fragment of 122 bp length. A 5' Rox labelled fluorescent anchor 
probe (acceptor) was designed with a phosphorylated 3' tail. The 
fluorescent sensor probe (donor) was designed with a 3' Fam 
modification (Table S3 in File SI). At the 3'-end of the sensor 
probe the donor fluorescent molecule is excited at a wavelength of 
533 nm. The fluorescent acceptor molecule at the 5'-end of the 
anchor probe receives the emitted energy from the donor. The 
emission of the fluorescence signal from the acceptor molecule is 
measured by the LightGycler 480 instrument using a filter 
combination of 483-610 nm. The sensor probe was designed to 
match the G-aUele of SNP ss83 18840 14. In the case of a 



mismatch due to presence of the G-aUele, a lower melting 
temperature is detected. The LightGycler 480 software was 
adjusted to automatically detect genotypes based on melting curve 
results. All genotypes were manually checked for genotyping 
errors, deciphered in triplicates via Sanger sequencing. Homozy- 
gous G/G individuals showed a fluorescence peak between 483 
and 610 nm wavelength at 58°G, whereas homozygous G/G 
individuals showed one at 64°C, and heterozygous fish showed 
two peaks (Figure SI in File SI). The melting curve analysis 
comprised of an initial deiiaturation step (95°C for 1 min), a step 
rapidly lowering the temperature to 40°G and holding for 30 sec, 
and a heating step slowly increasing the temperature up to 80°G 
under a continuous measurement of fluorescence ( 1 acquisition/ 
°G). In total 337 individuals were investigated using the FRET 
assay, in order to obtain a total of 430 genotypes at SNP 
ss831884014 (including the before mentioned Sanger-sequenced 
individuals). Briefly, the PGR was run in a 25 |J,1 reaction volume 
consisting of 1 xPGR buflFer containing MgGl2, 0.4 pmol of each 
primer and each FRET probe, 10 mM dNTPs, 2 U FastStart Taq 
DNA polymerase, HgO and 20 ng of DNA. All PGRs were carried 
out in 96 well plates in a LightGycler 480 with an initial 
denaturation at 95°G for 10 min, followed by 35 cycles of 92°G for 
30 s, 56°G for 30 s and 72°G for 30 s. 
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Figure 2. Gene structure, pairwise linicage disequilibrium [r^ between SNPs, and functional annotation of the Nile tilapia amh gene 
derived from the EnsembI database, a) Gene structure of the amh gene derived from Scaffold GL831 234.1. Exons 1-7 (El -7) are represented by 
filled grey boxes. Blue and red dashed lines show positions of polymorphic SNPs in the amh gene. SNP S5831884015 (C>T) is located -225 bp 
upstream of the start codon ATG, SNP ss831884018 was located in the intron between exon 6 and 7. Variants ss831884014 and ss831884019 were 
found in exons 6 and 7. SNP 5s831884014 (G>C) and S5831884019 (C>T) were missense mutations leading to amino acid changes from glutamic acid 
to glutamine (codons Gaa/Caa) and alanine to valine (codons gCg/gTg). Table 1 deals with the functional annotation of all six SNPs. b) Pairwise 
linkage disequilibrium heat map of four polymorphic allelic variants in the Nile tilapia amh gene, c) Partial DNA to protein translation of exon 6 in the 
Nile tilapia amh gene, depicting a non-synonymous amino acid substitution at position 376 of the putative amh-protein (indicated by black arrow) 
coded by an allelic variant 55831884014 of scaffold GL831 234.1. d) Partial comparative alignment of the amh genomic exon 6 sequences for six teleost 
species derived from the EnsembI database, e) Partial comparative alignment of the amh protein sequence for six teleost species derived from the 
EnsembI database. 

doi:1 0.1 371 /journal.pone.01 04795.g002 



Statistical analysis 

For the segregating SNPs (n = 4) the gene diversity, tlie allele- as 
well as genotype frequencies were analysed using the SAS/ 
Genetics 9.3 software (SAS Inst., Inc., Gary, NG, USA). Phased 
parental haplotypes were reconstructed using the software package 
PLINK tor 93 animals treated at high temperature [38]. The 93 
individuals comprised three families including the corresponding 
three dams and sires. /?^-values were calculated in SAS/Genetics 
9.3 software as a measure of linkage disequilibrium between 
segregating SNPs [39,40]. Associations of segregating SNPs with 



the phenotypic sex of the fish were investigated, fitting a 
generalized linear model (GLM) with binominal error distribution 
and logit function in SAS version 9.3. First, a stepwise logistic 
regression analysis was carried out to detect associations between 
any of the SNP genotypes coded as 0 (GC), 1 (GC), or 2 (GG) and 
sex coded as a binary trait (0 = male, 1 = female). Second, a gene 
substitution model was applied to detect sex specific effects and 
those of significant SNP driven from the stepwise regression. A 
stepwise logistic regression revealed a significant effect of SNP 
ss83 188401 4 and ss8 3 1884018 on the probability of developing 
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a male functional phenotype (stepwise logistic regression, binary 
logit, chi-square 9.81, p<0.0017). Both SNPs showed a significant 
effect of equal magnitude on the male proportion in progenies 
treated at high temperature as they were in full linkage 
disequilibrium {LD, r^= 1, no recombination event occurred 
between the two SNPs). However, SNP ss831884018 was not 
further considered as the potentially causal variant due to its 
position in the non-coding region (intron 6; Table 1). The 
following model was applied: 
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Logit 
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where Tirs is the probability of obtaining the male phenotype, (p is 
the overall mean effect, otr is the frxed effect of the treatment 
(levels: 28°C or 36°C); Y; is the linear regression coefficient for 
association between the probability of obtaining the male 
phenotype and SNP ss83 18840 14 genotype, X,,, is the effect of 
SNP SS831884014 genotype (levels: 0 (C/C), 1 (G/C), 2 (G/G)); 
Y2 is the linear regression coefficient for the fixed interaction effect. 
For derivation of the gene substitution effect only the significant 
parameters were considered in the final model. Since the inverse 
link is nonlinear and does not give an equal substitution effect, in 
addition the model parameters were estimated by a general linear 
model (GLM) using a normally distributed response variable with 
an identity link function. 

Additive and dominance effects of alleles were calculated 
according to the following model: 
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where ti^s is the probability of obtaining the male phenotype, <p is 
the overall mean effect, a,- is the fixed effect of the treatment 
(levels: 28°C or 36°C from 10-20 dpf); is the fixed effect of SNP 
SS831884014 genotype (levels: 0 (C/C), 1 (G/C), 2 (G/G)); 0(P„ is 
the fixed effect of interaction between genotype and treatment. 
First, least square means were estimated on the logit scale and then 

exp(x/?) 
[l+exp(x/J)l 

to the original scale (probability) applying the least square means 
(LSmeans) statement. Significant differences between least square 
means were tested using a t-test by inclusion of the PDIFF option 
in the LSmeans statement. Standard errors of least square means 
were calculated as described by [41]. Significant deviation of 
estimates of dominance effects from zero was tested using a t-test 
as described in [42]. 

Results 
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Phenotypes: Sex ratios in controls and groups treated at 
high temperature 

Temperature had a significant effect on the sex ratio of the 
genetically female (XX) line of Nile tilapia consisting of 1 2 families 
investigated in the present study. The overall male percentage 
when reared at control temperature of 28°C, was 2.3%. Losses 
from 10 dpf until sexing in control groups were 8.6%. In the 
corresponding fuU sib groups treated at high temperature, the 
male percentage was 40.4%, while mortalities were 4.4%. 

In the genotyped sample, which was derived from the before 
mentioned genetically female (XX) line of Nile tilapia consisting of 
12 families, a male proportion of 20% in controls, and 53% in the 
group treated at high temperature was observed. Similarly, the 
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temperature treatment significantly afiected tlie proportion of 
phenotypic males (p<0.0001). 

Functional annotation of amh variants 

In total, sequencing of the amh gene revealed six SNPs 
{ss831884014 to ss831884019, Table 1, Figure 2). SNP 
ss831884015 (OT) was detected in the 5' upstream region, 
located —225 bp upstream of the start codon, and was classified as 
intergenic modifier. Two of the identified SNPs were located in 
introns (intron 4 SNP ss831884016, and intron 6 SNP 
ss83 1884018). Finally, three variants were assigned to exons 5 
(ss83 1884017), 6 {ss83 1884014) and 7 {ss83 1884019). SNP 
ss83 1884017 was a synonymous mutation at position 696 of the 
coding sequence, corresponding to position 232 of the protein. 
Moreover, SNP ss831884014 (G>C) and s.s831884019 (OT) 
were missense mutations leading to amino acid changes from 
glutamine to glutamic acid (codons Gaa/Caa) and alanine to 
valine (codons gCg/gTg). 

Frequency distribution of geno- and haplotypes 

Four out of six variants (SNP ss831884015, ss831884014, 
ss831884018, and ss831884019) segregated in the investigated 
93 individuals from the genetically female Nile tilapia population 
(Table 2). Moreover, four phased haplotypes from parent to 
offspring were reconstructed (C-C-A-C, C-G-G-C, C-G-G-T, T- 
G-G-T). Variants ss831884014 and ss831884018 were in fuU 
linkage disequilibrium [LD). SNPs s.s831884015 and 
ss831884019 also exhibited a high degree of LD, showing an 
r^-value of 0.7. The LD between SNPs ss831884014 and 
SS831884019 as weU as between ss831884018 and 
ss83 1884019 was intermediate, showing an r~ -value of 0.54. In 
contrast, LD between SNPs ss831884015 and ss831884014, as 
well as ss831884015 and ss83 188401 8 was lower, showing an re- 
value 0.21 (Figure 2, Table S4 in File SI). 

Association of amh variant ss831884014 with autosomal 
and temperature-dependent sex reversal 

A fluorescence resonance energy transfer (FRET) based 
genotyping assay was developed for SNP S.s831884014 and 
validated on a representative sample of fish, comprised of 125 
individuals reared at control temperature (28°C) and 305 
individuals reared at 36°C from 10-20 dpf The statistical analysis 
revealed a highly significant effect of the treatment and genotype 
on the phenotypic sex (Table 3). In contrast, no significant 
genotype xtreatment interaction, and thus no difference in 
association between fish reared at control or high temperature, 
was detected (interaction term: df = 424, F=1.83, p = 0.1618. 



Table 3). Furthermore, using a reduced model including only the 
two main factors SNP and treatment as well as applying a 
Bonferroni correction resulted in significant differences between 
the group treated at high temperature and the control (p<0.0001), 
and SNP-genotypes (CC vs. GG; p<0.0007; GC vs. GG; p< 
0.0067), as no significant interaction between the SNP and 
treatment was observed (Table 3). 

Generally, the genotype CC showed a more pronounced effect 
on the percentage of males. Thus genotypic values over both 
treatments, given as LSmeans, were 0.4224 for the CC- compared 
to a significantly lower value of 0.1699 for the GG-genotype (p< 
0.0007). The genotypic LSmean of the heterozygous GC-genotype 
was 0.3604, which differed significantly from the genotype GG 
(p<0.0067). 

Although no significant interaction between treatments was 
detected, the genotypic LSmeans derived for both the group 
treated at high temperature and the control group are given 
separately in figure 3. The CC-genotype exhibited the highest and 
the GG-genotype the lowest genotypic LSmean male proportion, 
irrespective of the treatment. Thus, homozygous CC individuals 
showed an LSmean proportion of males (57.8%) when treated at 
high temperature. In comparison, homozygous CC specimens 
kept at 28°C exhibited LSmeans of 38.4%. The LSmean male 
proportions dropped to 31.9% and 5.0% for individuals possessing 
the homozygous genotype GG, in the group treated at high 
temperature and the control group, respectively. 

Gene substitution effect, homozygous additive allele 
effect, and dominance effect of amh variant ss831 884014 
for autosomal and temperature-dependent sex reversal 

An analysis of covariance was carried out to estimate the gene 
substitution effect, where the effect of the SNP-genotype was 
considered as the regression term (equation 2). The estimated gene 
substitution effect was 12% for both treatments. Genotypic values 
derived from the corresponding regression model were 63%, 50%, 
and 38% for genotypes CC, GC, and GG in fish treated at high 
temperature. In contrast to this, lower genotypic values of 32%, 
19%, and 7% were observed in fish kept at the control 
temperature of 28°C for genotypes CC, GC, and GG, respective- 
ly. Moreover, the genotype effect of SNP ss831884014 on the 
male proportion was significantly higher in groups treated at high 
temperature, yielding 31% higher male ratios (LSmeans) on 
average of the three genotypes (Figure 3). Furthermore, the 
homozygous additive effect of the C-allele was estimated, which 
deals with the increase in male percentage contributed by each 
copy of the C-aUele. This effect was calculated as the mean 
between the two homozygous genotypes (a = Hcc~liGG)/2). In 



Table 2. Allele frequencies, polymorphism information content (PIC), heterozygosity and allelic diversity for four allelic variants in 
the amh gene detected in a genetically all-female (XX) Nile tilapia population. 





SNP locus 


n individuals 


Allele frequency 


PIC 


Heterozygosity 


Allelic diversity 


ss83 1 88401 S 


93 


0.90 (C) 


0.1595 


0.1935 


0.1748 


ss83 188401 4 


93' 


0.71 (C) 


0.3245 


0.4194 


0.4075 




430^ 


0.58 (C) 


0.3677 


0.5186 


0.4856 


ss83 18840 18 


93 


0.71 (A) 


0.3245 


04194 


0.4075 


ss83 188401 9 


93 


0.84 (C) 


0.2340 


0.3226 


0.2706 



(' Initially 93 individuals were Sanger-sequenced for tine amh gene. Tlie genotyped sample was comprised of three families including the corresponding three dams and 
sires; ^ In total 337 individuals were investigated using the FRET-assay, in order to obtain a total of 430 genotypes at SNP ss831884014 {including the before mentioned 
Sanger-sequenced individuals)). 
doi:l 0.1 371/journal.pone.Ol 04795.t002 



PLOS ONE I www.plosone.org 



6 



August 2014 | Volume 9 | Issue 8 | el 04795 



Amh Polymorphisms and Sex Reversal in Nile Tilapia 



Q- 

o 

OJ 
Q. 
OJ 



CP 



I 

O 

E 
o 

«4- 
U 

vo 
ro 

i/i 
> 

o 

00 



OJ 

ro 

OJ 
Q. 

E 

OJ 



OJ 
Q. 
>. 

O 
c 

OJ 



s 

00 



u 



CM T— 

O 

O 



<— U-i -— 



o vo m CO -— 
^ «^ rsj 

'^r >— rsi CO -— 





OJ 
Q. 






ent 


?not 


NP ( 


ent 


atm 


Pg( 


LD 
fO 


atm 


OJ 


SN 




tre 



>- on 



4^ 4-" l/l 



-J -5 



+ 

+ 
+ 



controls the homozygous additive allele effect showed a value of 
16.7% (S.E. = 6.9), whereas a value of 12.9% (S.E. = 4.04) was 
observed in fish treated at high temperature (Table 4). In addition 
to the additive effects of alleles, dominance is of theoretical and 
practical importance (i.e. in breeding). A significant dominance 
deviation of 1 1 .4% (S.E. = 5.6) was estimated for the fish treated at 
high temperature. In contrast, no significant dominance deviation 
was detected in individuals reared at control temperature (— 
4.01%, S.E. = 8.6). 

Discussion 

THapias have the potential to become the most important 
aquaculture species group world-wide [43]. They serve as a means 
for poverty alleviation and concurrentiy constitute a world-wide 
commodity, being sold as frozen fillets to export markets in the 
US, Europe, and Asia [44]. More than 50% of seafood consumed 
comes from aquaculture; tilapia is a major mainstay of this 
production. 

Tilapia differ from other aquaculture species insofar as they can 
thrive under poor conditions and on diets that contain low shares 
of fish meal [2], but a major sustainability constraint - i.e. 
hormonal sex-reversal - is yet to be overcome. The present study 
reports novel allelic variants in the Nile tilapia amh gene. One of 
the variants is associated with the sex and might enable marker- 
assisted selection. A strategy is proposed to increase the male 
proportion by the development of a straightforward genotyping 
assay combined with temperature treatment of fry. 

Association of amh variant ss831884014 with autosomal 
and temperature-dependent sex reversal 

The largest effect on the phenotypic sex of Nile tilapia was 
observed in homozygous carriers of the C-aUele at SNP 
ss831884014, a missense variant that leads to an amino acid 
change from glutamine to glutamic acid at position 376 of the 
putative mnh protein. Moreover, a ten-day temperature treatment 
of the sexually undifferentiated fry carrying the homozygous CC 
genotype significandy increased the probability of developing the 
male phenotype. Although the homozygous effect of the C-aUele 
was higher in the control groups (16.7% compared to 12.9% 
males), the average genotypic values were 31% higher (LSmean 
male ratios) in groups treated at high temperature. 

Two factors were identified which suggest that the C-aUele (in 
homozygous and heterozygous aUelic state) in our population has 
the potential to significantly increase the proportion of males: 1) 
Genotypic values of 63% and 32% males for the CC-genotype and 
38% and 7% for the alternative GG-genotype in the group treated 
at high temperature and control group, respectively; 2) A 
significant dominance deviation (+11.4% of males) of the 
heterozygous (CG) genotype in the fish treated at high temper- 
ature. 

Comparative genomics and functional role of the amh 
variant ss831884014 

A comparative analysis of the amh protein revealed that the G- 
aUele might be the ancestral aUele at locus ss831884014 (Figure 2, 
panel e). Species such as tilapia, stickleback, and medaka have the 
amino acid glutamic acid (E: coded by Gaa, Gag) in their protein 
sequence, whereas puffer fish and platy fish exhibit the acidic 
amino acid aspartate (D: coded by Gac, Gat). Hence, this might 
indicate that glutamine (Q), which is coded by the nucleotide 
codon Caa, is more recent in the Nile tUapia line investigated here 
and might have taken over a critical role in autosomal and 
temperature-dependent sex reversal. 
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Figure 3. Relationship between male proportion and genotype at amh variant ss831884014 in Nile tilapia, reared at 28°C or 36°C 
from 10 to 20 days post fertilisation. Linear least-square regression of the genotypic values (LSmean male proportion) for genotypes C/C, G/C, 
and G/G at locus S5831884014 of the amh gene (Scaffold GL831234.1). The genotypic values were derived from the regression models by either 
excluding the interaction between SNP and treatment (LSmeans 1: red and blue line) or by including the interaction between SNP and treatment 
(LSmeans 2: red and blue triangles) as effect class. Numbers of individuals per genotype were 26/79/20 and 1 14/144/47 for genotypes C/C, G/C, and 
G/G, for the control group (28"C) and the group reared at high temperature (36 C), respectively. 
doi:1 0.1 371/journal.pone.01 04795.g003 



The mnh gene was identified as a potentially sex-modifying cue 
in a number of earlier studies [15,16,18,22,32-34]. On the one 
hand, expression analysis confirmed amh as a functional marker of 
maleness at a precocious age of 3 dpf and later during sex 
differentiation [15,32-34]. On the other hand, evidence for sex- 
specific differences in the amh gene at the genomic level was 
lacking so far. An earlier identified variant [AM232733), similarly 
located in exon 6 of the gene, proved not to be a trigger for sexual 
fate in Nile tilapia [18]. Moreover, the Nile tilapia population 
investigated here (Lake Manzala population), showed no segrega- 
tion of variant ^M232 753. However, among 10 genes putatively 
involved in sex determination or differentiation, amh is located in 
the centre of an SD QTL on LG23 [15]. 

The role of the amh variant ss831 884014 in the 
continuum of major genetic, minor genetic, and 
temperature-dependent factors 

The present study does confirm - for the first time - that an 
allelic variant in the amh gene might be a major QTL for 
autosomal and temperature-dependent sex reversal in Nile tilapia. 
In the Patagonian pejerrey a functional duplication of the amh 
gene suggests that amhy may be the master sex-determining gene 
[45] , but our results for tilapia point in another direction: multiple 
interacting loci each partially contributing to the formation of the 
sexual phenotype. Although the allelic variant reported here exerts 
a large effect on the formation of the sexual phenotype under both 
rearing environments, not all sex phenotypes can be explained 
using a simple one-locus model. Recently, restriction site 
associated sequencing revealed the main sex-determining locus 



flanked by two SNPs [Oni23063 and Oni28137) and located on 
linkage group 1 [46] . The authors showed that in one family, sex- 
reversed males which were treated at high temperature carried the 
female genotype at SNP Oni23063 and Oni28137, but did not 
show an association with the temperature-dependent sex [46] . The 
underlying gene of the main sex-determining locus on LGl is 
unknown so far. Further studies are needed to verify if SNP 
ss83 1884014 on LG23 acts as autosomal sex modifying cue in this 
Nile tilapia line, which in addition or alternatively to the major 
genetic factor on LG 1 , described in [46] , increases the probability 
of developing the male phenotype during temperature-dependent 
sex reversal. The missing significant interaction between the SNP- 
genotype and treatment effect detected in the present study 
(Table 3) suggests that the same genetic network (i.e. the TGF-fi 
pathway) might be at least partially active during autosomal and 
temperature-dependent sex reversal in this Nile tilapia line, yet 
with a drastically pronounced effect under elevated rearing 
temperatures from 10-20 dpf 

Further research should now focus on the functional role of the 
identified allelic amh variant during autosomal and temperature- 
dependent sex reversal. The underlying gene networks, such as the 
TGF-P pathway, need to be investigated in greater depth, and 
their transcriptional, translational, and phenotypic consequences 
have to be evaluated to fuUy characterize the effect on the 
phenotypic sex described here. Targeted techniques such as RNA- 
seq [47], or morphoKno-mediated knockdown of known target 
genes might reveal the functional role of genes and variants [48] . 
Another avenue of enquiry focuses upon the specific biochemical 
properties of the amh hormone. The combined use of temperature 
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treatment and marker-assisted selection to introgress the C-allele 
at W? ss83 18840 14 into otiier strains of cultivated tilapia might 
help to increase the male proportion and therefore the overall 
productivity. Ultimately, increasing the yield through higher male 
proportions without counteracting reproductive traits would 
greatly contribute to increase the efficiency of Nile tilapia culture 
and protein availability in predominant small-scale aquaculture, 
while also minimizing negative environmental effects due to 
hormonal sex-reversal from expanding intensive aquaculture 
systems. 

Conclusions 

This study shows that an allelic variant in the aynh gene is a 
major QTL for autosomal and temperature-dependent sex 
reversal in a selected line of Nile tilapia. If the applicability of 
this marker in other populations proofs successful, a combined 
strategy of marker-assisted selection and temperature treatment 
might be beneficial to increase the male proportion in aquaculture 
stocks of O.niloticus. 

Supporting Information 

File SI Contains supporting Figures and Tables. Figure 
S 1 , Melt curve analysis output from the Lightcycler 480 system for 
genotypes at amh variant ss83 18840 14. Different genotypes are 
illustrated using dififerent colours: C/C = green curve, C/G and 
G/C = blue curve, and G/G = red curve. Homozygeous C/C 
individuals show a fluorescence peak between 483 and 610 nm 
wave length at 58°C, whereas homozygeous G/G individuals 
showed one at 64°C, and heterozygous fish showed two peaks. 
Table SI, Pedigree and sex ratios of the genetically female 
population reared at control (28°C) and elevated temperature 
(36°C) from 10 to 20 dpf Table S2, Forward and reverse primers 
tailed with a universal Ml 3 forward or reverse primer for 
bidirectional sequencing the amh gene in Nile tilapia. Table S3, 
Fret-primer for allelic variant 1690582 in the NUe tilapia amh 
gene, anchor and sensor probe sequences and positions on scaffold 
GL831234.1. Table S4, R^-measure of linkage disequilibrium 
between four segregating allelic variants in the amh gene of Nile 
Tilapia. The estimates were derived from a sample 93 tempera- 
ture-treated Nile tilapia individuals. Table S5, Raw data for the 
genetically female study population reared at control (28°C) and 
elevated temperature (36°C) from 10 to 20 dpf Table S6, 
Genotypes of four segregating SNPs in the atnh gene of 93 
individuals derived from three Nile tilapia families. 
(ZIP) 
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